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Abstract 

The  density  and  viscosity  of  propane  mixed  with  66/34  trans/cis- decahydronaphthalene  were 
measured  over  a  wide  range  of  temperatures  (323  -  423K),  pressures  (2.5-208  bar),  and 
compositions  (0-65  mol%  propane).  For  conditions  giving  two  phases,  the  composition  of  the 
dense  phase  was  measured  in  addition  to  the  density  and  viscosity.  The  modified  Sanchez- 
Lacombe  equation  of  state  (MSLEOS)  was  used  with  a  single  linearly  temperature-dependent 
pseudo-binary  interaction  parameter  to  correlate  the  phase  compositions  and  densities.  The 
compositions  and  densities  of  the  mixtures  were  captured  well  with  absolute  average  deviations 
between  the  model  and  the  data  of  5.3%  and  2.3%,  respectively.  The  mixture  viscosities  were 
computed  from  a  free  volume  model  (FVM)  by  using  a  single  constant  binary  interaction 
parameter.  Density  predictions  from  the  MSLEOS  were  used  as  input  mixture  density  values 
required  for  the  FVM.  The  FVM  was  found  to  correlate  well  with  the  mixture  viscosity  data 
with  an  absolute  average  deviation  between  the  model  and  the  data  of  5.7%. 
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1.  Introduction 

Organic  solvents  expanded  with  supercritical  fluids  (ScF)  have  been  investigated  as 
alternative  chemical  process  media  for  more  than  two  decades.  ScF  expanded  liquids  can 
provide  an  otherwise  difficult-to-achieve  range  of  solute  solubilities  obtained  by  mixing  a  good 
solvent  (the  organic  liquid)  with  a  poor  solvent  (the  ScF).  Recently,  ScF/solvent  combinations 
have  been  recognized  as  unique  and,  in  some  cases,  beneficial  chemical  reaction  solvents  [1]. 
For  mass  transfer-controlled  reactions,  the  addition  of  a  ScF  can  improve  reaction  rates  by 
increasing  mass  transfer  rates  in  the  fluid  phase  and  to  the  internal  surfaces  of  porous  catalysts. 
As  examples,  solvents  expanded  by  supercritical  and  subcritical  ScFs  have  been  used  in 
homogeneous  catalytic  oxidations  [2],  heterogeneous  catalytic  oxidations  [3],  and  heterogeneous 
hydrogenations  {4-7}. 

The  addition  of  a  ScF  to  polymer/organic  solvent  solutions  has  been  investigated 
extensively  as  a  technique  to  decrease  the  lower  critical  solution  temperature  (LCST)  allowing  a 
less  aggressively  attainable  phase  split  for  polymer  separation/purification  [8-14],  Recently,  the 
hydrogenation  of  polystyrene  (PS)  has  been  investigated  in  decahydronaphthalene  (DHN) 
expanded  by  supercritical  carbon  dioxide  (SCCO2)  [4,  5,  7],  Although  the  addition  of  scCC>2 
improved  the  hydrogenation  rate  under  many  circumstances,  unwanted  carbon  monoxide  (CO) 
was  produced  via  the  reverse  water  gas  shift  reaction  (RWGSR).  When  the  concurrent 


production  of  CO  is  unavoidable,  methanation  metals  like  Ni  or  Ru  can  be  used  to  mitigate  CO 
poisoning  of  the  hydrogenation  catalyst  [4,  7]. 

Alternatively,  propane  (C3)  can,  in  principle,  replace  CO2  in  this  reaction  to  eliminate  the 
formation  of  CO  entirely.  The  results  presented  here  are  germane  to  that  alternative  for  PS 
hydrogenation.  Due  to  the  high  cost  of  the  pure  isomers  of  DHN,  a  commercial  mixture  of 
trans/cis- DHN  produced  by  DuPont  (Decalin®)  was  used  for  our  studies.  The  composition  of 
commercially  available  DHN  samples  ranges  from  almost  100%  tnmv-DHN  to  an  isomeric 
mixture  with  as  little  as  40%  trans-DHN  [15].  Polystyrene  phase  separation  behavior  in  DHN 
has  been  reported  and  has  been  interpreted  in  light  of  the  properties  of  the  solvent  [16].  For  C3 
expanded  DHN,  fundamental  thermophysical  data  are  lacking  in  the  literature  and 
thermophysical  data  for  PS  in  C3/DHN  mixtures  have  not  been  published.  In  this  paper,  we 
report  thermophysical  data  for  mixtures  of  DHN  and  C3.  We  use  the  modified  Sanchez-Lacombe 
Equation  of  State  (MSLEOS)  to  fit  phase  density  and  composition  data  and  a  free  volume  model 
(FVM)  to  fit  the  mixture  viscosity  data.  This  data  is  of  significant  importance  in  analyzing  the 
performance  of  PS  hydrogenation  reactions  in  high-pressure  C3/DHN  mixtures. 

2.  Materials  and  Methods 

Research  grade  C3  (99.99%,  Airgas  National  Welders)  and  high  purity  66/34  wt% 
trans/cis- DHN  (99+%,  Sigma-Aldrich)  were  used  for  all  measurements.  All  chemicals  were 
used  as-received.  The  phase  behavior  and  densities  of  C3/DHN  mixtures  were  measured  in  a 
variable-volume  view  cell  housed  within  a  thermostatted  (±  0.2  °C)  oven.  A  schematic  diagram 
of  the  experimental  apparatus  is  shown  in  Figure  1.  Pressure  was  transmitted  from  the  sample 
fluid  through  a  chemical  seal  (L990.34,  WIKA  Instrument  Corp.)  to  a  pressure  transducer 
(PX309-10KGV,  Omegadyne,  ±  1.7  bar)  mounted  outside  of  the  oven.  Silicon  oil 


(Phenylmethyl  Siloxane,  Dow  Coming)  was  used  as  the  pressure  transmitting  fluid.  Linear 
displacement  of  the  piston  was  measured  using  a  micrometer  (CD-6”CSX,  Mitutoyo).  A 
calibration  curve  of  volume  versus  displacement  was  made  by  adding  a  known  amount  of  CO2 
and  measuring  the  displacement  of  the  piston  at  temperatures  and  pressures  where  the  densities 
were  known.  Pure  component  and  mixture  densities  were  measured  using  the  calibration  curve 
along  with  the  mass  of  the  sample  loaded  into  the  view  cell.  All  wetted  surfaces  were  Teflon™ 
or  316  stainless  steel. 


Figure  1  Schematic  of  phase  equilibria  apparatus  [TI:  temperature  indicator,  PI:  pressure 
indicator,  TIC:  temperature  indicator/controller,  DG:  displacement  gauge] 

Viscosities  were  measured  using  a  SPL440  high-pressure  viscometer  (Cambridge 
Applied  Systems,  Cambridge,  MA)  equipped  with  a  variable-volume  view  cell  and  an  external 
recirculation  pump  (Micropump  GAH.X21,  Idex,  Oak  Harbor,  WA)  used  for  mixing.  A 
schematic  of  the  experimental  apparatus  is  shown  in  Figure  2.  The  external  gear  pump  head  was 
rapped  with  insulation  to  minimize  heat  losses  during  recirculation.  The  pressure  was 
maintained  well  above  a  predetermined  pressure  to  ensure  a  single  fluid  phase.  The  mixture  was 
assumed  homogeneous  when  convergence  of  the  measured  viscosity  was  obtained  after  checking 


the  viscosity  at  different  intervals  during  recirculation.  Once  the  mixture  was  assumed 
homogeneous  the  2- way  and  3 -way  valves  were  closed  to  the  recirculation  loop  to  isolate  the 
measurement  fluid  in  the  oven  from  the  fluid  outside  the  oven  in  the  recirculation  pump. 
Precautions  were  taken  to  ensure  phase  separation  did  not  occur  between  measurements  at  a 
given  composition.  The  viscosity  of  the  fluid  is  measured  by  recording  the  time  it  takes  the 
piston  to  move  between  the  two  magnetic  coils.  This  time  is  converted  to  a  value  of  viscosity  in 
centipoise  (cP)  by  an  internal  calibration  supplied  by  the  manufacturer.  Before  each  viscometer 
experiment  the  viscometer  apparatus  was  leak  checked  by  charging  with  carbon  dioxide  and 
pressurizing  to  346  bar.  The  viscometer  was  evacuated  and  charged  with  propane  three  times  to 
remove  residual  carbon  dioxide.  For  an  actual  viscosity  experiment,  C3  was  charged  to  the  cell 
with  a  high  pressure  syringe  pump  (500D,  Teledyne  Isco).  The  mass  of  C3  added  to  the  view 
cell  was  computed  from  the  change  in  volume  of  the  syringe  pump  at  the  temperature  and 
pressure  of  the  syringe  pump  using  the  known  densities  of  C3.  A  sample  bomb  was  used  to 
measure  the  accuracy  of  the  C3  addition  procedure.  The  difference  between  the  calculated  mass 
from  the  syringe  pump  and  gravimetrically  measured  mass  added  to  the  sample  bomb  was  at 
most  0.1%.  The  estimated  standard  uncertainties  for  temperature,  pressure,  and  C3  mole  fraction 
were  0.05  K,  0.06  bar,  and  0.005  mole  C3/mol,  respectively.  The  estimated  combined  expanded 
uncertainty  for  viscosity  and  density  were  0.027  cP,  and  0.002  g/cm3,  respectively. 


viscometer 


Figure  2  Schematic  of  high  pressure  viscometer  experimental  apparatus 
The  saturated  liquid  density  and  vapor  pressure  of  C3  are  well  established  [17,  18]. 
Therefore,  we  used  C3  to  validate  the  operation  of  the  view  cell  for  vapor  pressure  and  density 
measurements.  The  measured  vapor  pressures  are  compared  to  literature  values  in  Figure  3.  The 
vapor  pressure  was  measured  by  recording  the  pressure  when  the  first  bubble  appeared  during  a 
decompression  experiment.  The  process  was  repeated  several  times  to  ensure  reproducibility. 
The  vapor  pressure  was  also  measured  by  recording  the  pressure  when  the  cell  was 
approximately  50%  liquid  and  50%  vapor  and  was  allowed  to  equilibrate  for  approximately  one 
hour.  There  was  negligible  difference  between  the  two  methods  so  we  have  not  separated  values 
accordingly  in  Figure  3.  Volumes  were  measured  by  precision  micrometry  as  described 
previously  in  the  text.  Saturated  liquid  densities  of  C3  were  measured  at  the  bubble  point  and  are 
plotted  in  Figure  4. 


40 

35 

30 

25 

20 

15 

10 

5 

0 

f 

re  3 


i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

o  Haynes  et  al.  • 

a  Younglove  et  al.  ° 

•  This  work 


o 


o 


•° 


o 


o 


& 


o 


100  150  200  250 

Temperature  (K) 


300 


350 


400 


Comparison  of  measured  and  reported  C3  saturated  vapor  pressure  data 


E 


o 

3 

>* 

c 

a> 

Q 


0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 

50  100  150  200  250  300  350  400 


o 


o 


o 


o 


o« 


o 


o 


•a 


°A 


o  Haynes  et  al. 

•  Youngblood  et  al. 
a  This  work 

_l _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ l_ 


I  . 


Temperature  (K) 


Figure  4  Comparison  of  measured  and  reported  C3  saturated  liquid  density  data 

3.  Phase  equilibria  and  density  model 

The  lattice-fluid  model  proposed  by  Sanchez  and  Lacombe  [19,  20,  21],  expresses  the 
reduced  density  as  a  function  of  the  reduced  temperature  and  the  reduced  pressure  (Equation  1). 
Their  EOS  has  been  used  to  describe  the  thermodynamic  properties  of  many  complex  systems, 
including  organic  liquid/ScF  mixtures  [22],  polymer/organic  liquid  solutions  [16,  23],  and 
polymer/organic  liquid/ScF  solutions  [14,  24],  The  Sanchez-Lacombe  equation-of-state 
(SLEOS)  uses  three  pure  component  parameters  to  characterize  the  phase  behavior.  Typically, 
the  parameters  are  obtained  by  fitting  the  EOS  to  pure  component  saturated  vapor  pressure  data 


[20], 


P  =  -T  ln(l-p)  + 
The  reduced  properties  are  defined  as: 
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Neau  showed  that  the  SLEOS  can  give  thermodynamically  inconsistent  phase  equilibria 
predictions.  This  occurs  when  the  SLEOS  is  combined  with  certain  commonly  used  mixing 
rules.  In  particular,  inconsistent  results  are  obtained  when  mixing  rules  are  computed  using 
literature  chemical  potentials  derived  on  the  basis  of  the  Helmholtz  energy  obtained  from  the 
partition  function  of  mixtures  [22],  More  specifically,  the  reference  chemical  potential  at  the 
ideal  gas  state  leads  to  compositional  dependence  for  commonly  used  mixture  rules  and  hence  is 
thermodynamically  inconsistent.  Neau  suggested  instead  the  use  of  fugacity  coefficients  derived 
from  the  SLEOS.  Her  approach  gives  thermodynamically  consistent  phase  equilibria  predictions 
and  her  results  correlate  better  to  binary  data  for  C02/n-hexane  than  did  predictions  using 
literature  chemical  potentials  with  the  same  mixing  rules. 

Following  Neau,  fugacity  coefficients  are  calculated  from  the  expression  [23], 
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The  corresponding  compressibility  factor  z  is  calculated  from  the  EOS  via  the  equation, 
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The  commonly  used  “  ktj ,  /y  ”  mixing  rules  given  by  McHugh  and  Krukonis  [25]  were  used  in  our 
work  and  are  defined  as  follows: 
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The  partial  derivatives  of  o*  and  s*  for  the  McHugh/Krukonis  mixture  rules  are  computed  from 
Equations  8  and  9,  respectively  [25], 


The  SLEOS  pure  component  parameter  values  obtained  from  the  literature  for  C3,  trans- 
DHN,  and  czs-DHN  are  shown  in  Table  1.  For  computational  simplicity,  a  pseudo-compound 
comprising  66/34  wt%  tra«,s'/c/.v-DHN  was  used.  The  pseudo-compound  pure  component 
parameters  were  obtained  by  application  of  the  mixing  rules  to  the  two  isomers  a  priori  in  the 
phase  equilibria  calculations.  In  doing  this  the  binary  interaction  parameters  for  the  two 
components  were  assumed  to  be  zero.  Ternary  phase  equilibria  calculations  for  this  system 
entail  little  added  complexity.  However,  by  treating  the  ternary  system  (C3,  ?ra«.v-DHN,  and  cis- 
DHN)  as  a  pseudo-binary  system  (C3  and  DHN),  a  single  gas/solvent  binary  interaction 
parameter  can  be  used  for  a  gas/solvent/solvent  calculation.  If  the  EOS  accurately  describes  the 


pseudo-binary  gas/solvent  PVT  behavior,  accurate  predictions  for  the  polymer  solutions  are 
possible  by  employing  the  pseudo-binary  simplification. 

Table  1  Reported  SLEOS  pure  component  parameters  [12]. 


Compound 

e*  (bar  cm3 /mol ) 

v*  ( cm3/mol ) 

r 

C3 

30836.1 

9.84 

6.50 

rr<ms-DHN 

51630.6 

16.37 

9.03 

cA-DHN 

52509.2 

15.74 

9.15 

66/34  wt%  Trans/c is -DHNa 

51927.2 

16.15 

9.07 

a  Calculated  from  pure  component  values  using  mixing  rules 

4.  Viscosity  model 

Tests  of  the  free-volume  viscosity  model  (FVM)  developed  by  Allal  et  al.  [26,  27], 
derived  from  a  theoretical  and  physical  background,  have  been  reported  extensively  in  the 
literature.  Studies  have  been  reported  for  both  pure  fluids  [28-32]  and  for  mixtures  [33-37].  The 
FVM  is  especially  useful  for  describing  the  behavior  of  dense  fluids  (>200  kg/m3).  The  model 
uses  the  three  pure  component  parameters  ( i ,  a  and  B)  to  correlate  the  pure  fluid  viscosity.  The 
parameter  a  is  characteristic  of  the  barrier  energy  a  molecule  must  exceed  in  order  for  self- 
diffusion  to  occur.  The  characteristic  length  parameter  t  accounts  for  the  size  of  the  molecule. 
The  parameter  B  characterizes  the  free  volume  overlap  where  the  free  volume  in  a  pure  state  is 
defined  as  the  difference  between  the  molar  volume  and  the  hard-core  volume  of  the  molecule, 
and  arises  when  the  Doolittle  theory  is  incorporated  in  the  derivation  of  the  FVM  [27],  With  the 
incorporation  of  an  appropriate  dilute  gas  viscosity  model,  the  FVM  can  be  extended 
successfully  to  predict  the  viscosities  of  gases,  fluids,  and  their  mixtures  when  densities  are  less 
than  200  kg/m3  [27],  The  Allal  et  al.  expression  for  the  viscosity  takes  the  form  [27,  28], 


The  binary  interaction  parameter  for  viscosity,  kay,  was  added  to  the  mixture  rules  to 
improve  the  fit  of  the  FVM  to  experimental  data.  The  binary  interaction  parameter  for  viscosity 
helps  to  mitigate  the  detrimental  impacts  on  curve  fitting  of  inaccuracies  in  the  fluid  densities 
needed  as  inputs  to  the  viscosity  model,  inaccuracies  in  pure  component  viscosity  predictions, 
and/or  any  shortcomings  of  the  mixing  rules  used  to  account  for  the  binary  interactions  of  the 
fluids. 


The  FYM  pure  component  parameters  for  C3,  //vms-DHN  and  c/.s-DHN  obtained  from 
the  literature  are  shown  in  Table  2.  Values  of  the  pure  component  parameters  for  a  variety  of 
other  hydrocarbons,  alcohols,  refrigerants,  carbon  dioxide,  and  water  can  be  found  elsewhere 
[26].  Parameter  values  for  alkoxyethanols  [30],  polyethers  [31],  and  polyalkylene  glycol 
dimethylethers  have  also  been  reported  [29,  38], 

Table  2  Reported  FVM  pure  component  parameters 


Compound 

a  (j  mol  1m3kg  1) 

f  (nm) 

B  X  103 

C3  [28] 

59.4963 

0.0847825 

7.392 

trans-DHN  [29] 

100.7 

0.09087 

10.51 

czs-DHN  [29] 

103.8 

0.07493 

12.03 

5.  Results  and  discussion 

5.1  Density  and  phase  equilibria 

Tobaly  et  al.  reported  spectroscopically  measured  phase  composition  data  for  C3 /cis- 
DHN  mixtures  [38],  We  established  phase  compositions  at  392  K  for  C3  in  civ-DHN  and  for  C3 
in  the  66/34  wt%  trans/cis- DHN  for  comparison  to  their  data.  The  good  agreement  between  our 
view  cell  measurements  and  those  by  Tobaly  and  coworkers  is  illustrated  in  Figure  5.  The 
solubility  of  C3  is  slightly  higher  in  tnmv-DHN  than  in  c/v-DHN.  The  trans  isomer  has  a  lower 
saturated  liquid  density  when  compared  to  the  cis  isomer.  Consequently,  tnmv-DHN  is  expected 
to  have  a  higher  free  volume  than  does  cxs-DHN  with  concomitantly  higher  C3  solubility.  In 
addition,  the  conformational  difference  between  the  isomers  could  impact  C3  solubility  but  there 
is  no  rationale  for  choosing  higher  C3  solubility  in  one  of  the  isomers  based  on  conformational 
preferences.  Similar  behavior  has  been  observed  with  the  solubility  for  hydrogen  in  the  two 
isomers  of  DHN  [39]  and  we  have  also  observed  similar  differences  with  CO2  solubility  in  DHN 
[40].  Nevertheless,  based  on  the  experimental  findings,  the  solubility  of  C3  is  expected  to  be 
slightly  higher  in  the  isomeric  mixture  than  in  pure  c/'s-DHN. 


Mole  Fraction  (mol  C3/mol) 

Figure  5  Comparison  of  measured  C3  solubility  in  cz's-DHN  and  66/34  wt%  trans/cis- DHN 
and  reported  results  for  C3  in  cis-DHN  at  392  K 

Bubble  points  for  mixtures  of  C3  (5-65  mol  %)  in  66/34  wt%  trans/cis-DHN  were 
measured  at  323,  373,  and  423  K.  The  SLEOS  was  fit  to  the  bubble  points  using  a  linearly 
temperature  dependent  binary  interaction  parameter,  kc3-DHN,  in  Equation  6.  The  value  of  l..  was 

set  to  zero.  The  P-x  plot  comparing  the  experimental  data  and  the  curves  computed  using  the 
SLEOS  is  shown  in  Figure  6,  which  also  includes  the  formula  for  the  temperature  dependent 
binary  interaction  parameter  kc3-DHN-  The  curves,  calculated  by  using  the  pseudo-binary  mixture 
approach,  clearly  give  an  acceptable  fit  to  the  data  with  an  absolute  average  deviation  (AAD)  of 
5.3%.  The  experimental  uncertainty  for  the  bubble  point  measurements  was  3.5%.  We  also 


fitted  the  data  using  a  three-component  model.  Since  deconvoluting  the  civ-DHN/C3,  trans- 
DHN/C3,  and  cis- D  H N/tnmv-  D  H  N  pair-wise  interactions  was  problematic,  we  assumed  the 
C3/DHN  binary  interaction  parameters  were  equal  (kc3-c/mHN  =  kc3-?ra«iDHN)  for  the  three- 
component  calculation.  No  significant  difference  between  the  three-component  and  two- 
component  pseudo-binary  fit  to  the  bubble  points  was  found. 

The  single  binary  interaction  parameter  kc3-DHN  was  calculated  for  the  pseudo-binary 
mixture  of  C3  in  a  pseudo-component  comprising  66/34  wt%  trans/cis- DHN.  While  this 
parameter  is  specific  to  the  66/34  wt%  mixture  of  trans/cis  isomers  and  does  not  necessarily 
apply  to  other  pseudo-component  ratios,  our  ability  to  fit  P-x  for  both  c/s-DHN  and  the  66/34 
pseudo-binary  mixture  suggests  that  the  value  of  kc3-DHN  obtained  can  be  used  with  good  results 
for  other  trans/cis  combinations. 


Mole  Fraction  (mol  C3/mol) 

Figure  6  C3  solubility  in  pseudo  DHN  (66/34  wt%  trans/cis- DHN)  from  this  work;  solid 
lines  pseudo-binary  SLEOS  results 

Table  3  Tabulated  data  for  C3  solubility  in  66/34  wt%  trans/cis- DHN 

^  (mol  C3/mol)  PM  (bar),  323  K_ PM  (bar),  373  K_ P„  (bar),  423  K 


0.05 

2 

3 

3 

0.15 

3 

7 

10 

0.25 

5 

11 

18 

0.35 

7 

15 

25 

0.45 

9 

19 

32 

0.55 

10 

23 

40 

0.65 

12 

27 

48 

Miyake  and  coworkers  investigated  the  influence  of  stereoisomers  of  DHN  on  PVT 
properties  using  a  vibrating- tube  densimeter  [41].  Miyake  and  coworkers  tested  three  density 


models  against  their  data  for  pure  tnmv-DHN,  ch-DHN,  and  their  mixtures.  The  Peng-Robinson 
equation-of-state  (PREOS),  the  perturbed  chain-statistical  associating  fluid  theory  (PC-SAFT), 
and  the  molecular  dynamics  (MD)  approach  by  Miyake  et  al.  yielded  similar  absolute  average 
deviations  of  less  than  1%  for  294  points  at  7  different  compositions,  including  the  two  neat 
isomers,  three  temperatures  of  303.15,  323,  and  343.15  K,  and  14  different  isobars  over  the  range 
l-to-650  bars. 

We  have  compared  SLEOS  predictions  to  the  entire  data  set  of  Miyake  et  al.  and  their 
data  and  the  predictions  at  1  bar  and  650  bars  in  Figure  7.  For  the  set  of  294  data  points,  the 
AAD  obtained  using  the  SLEOS  was  0.38%.  The  SLEOS  compressed  liquid  density  predictions 
for  the  pure  isomers  mms-DHN  and  c/s-DHN  and  the  trans/cis-DHN  mixtures  capture  the  data 
qualitatively,  but  the  biased  error  (0.24%)  is  larger  than  the  experimental  error  (—0.06%)  reported 
by  the  authors.  Even  though  the  error  between  the  model  and  their  data  is  on  the  order  of  ten 
times  more  than  the  authors’  experimental  error,  the  model  predictions  compared  to  their  data  are 
still  more  than  adequate  for  engineering  design  work.  Indeed,  the  SLEOS  predictions  to  their 
data  are  better  than  the  three  models  tested  by  the  authors. 


Mole  Fraction  (mol  c/s-DHN/mol) 

Figure  7  SLEOS  density  calculations  from  this  work  for  DHN  mixtures,  data  from  Miyake 
et  al.  [41]  (unfilled  triangles- 1  bar,  343.15  K,  unfilled  squares- 1  bar,  323  K,  unfilled  circles- 1 
bar,  303.15  K,  filled  triangles-600  bar,  343. 15K,  filled  squares-600  bar,  323  K,  filled  circles-600 
bar,  303.15  K) 

The  saturated  liquid  densities  of  the  pseudo-binary  mixtures  of  C3  in  66/34  wt% 
trans/cis- DHN  were  measured  at  the  bubble  point  of  each  mixture  at  323K,  373  K,  and  423  K. 
In  Figure  8,  densities  computed  using  the  SLEOS  are  compared  to  the  experimental  solution 
densities.  There  does  appear  to  be  some  systemic  bias  for  the  density  predictions  but  the  AAD 
was  0.92%  and  the  experimental  uncertainty  was  estimated  as  0.83%. 
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Figure  8  Saturated  liquid  density  of  C3  in  66/34  wt%  trans/cis- DHN  mixtures  vs.  C3 
concentration  and  SLEOS  calculation  results  from  this  work 

Table  4  Tabulated  data  for  saturated  liquid  density  of  C3  in  66/34  wt%  trans/cis- DHN  mixtures 


Xp,  (mol  C3/mol) 

0.05 

0.15 

0.25 

0.35 

0.45 

0.55 

0.65 


P*  (g/mL),  323  K 

0.845 

0.832 

0.816 

0.787 

0.758 

0.728 

0.671 


p^lg/mL),  373  K 

0.813 

0.791 

0.759 

0.741 

0.708 

0.678 

0.638 


P*  (g/mL),  423  K 

0.769 

0.741 

0.714 

0.692 

0.661 

0.611 

0.565 


To  better  understand  the  volume  expansion 


(VE)  behavior  of  C3/DHN  mixtures, 


Equation  14  is  used  to  calculate  the  VE,  defined  as  the  specific  volume  of  the  C3/DHN  mixture, 


VMIX,  at  a  given  temperature  and  pressure  divided  by  the  specific  volume  of  the  DHN  mixture, 
VDHN,  at  the  same  temperature  and  pressure, 


VE  = 


V 


MIX 


DHN 


(T,P) 


(14) 


The  results  for  the  VE  of  the  saturated  liquid  and  the  compressed  fluid  at  one  isotherm  are  shown 
in  Figure  9. 

Two  other  VE-P  isotherms  were  measured  (323  K  and  423  K)  and  similar  behavior  was 
observed.  The  maximum  VE  measured  of  1.385  occurs  at  423  K  and  at  the  bubble  point  for  a 
mixture  of  65%  C3  in  DHN.  The  SLEOS  predictions  for  the  single-phase  mixture  densities 
capture  the  qualitative  behavior  of  the  data  reasonably  well  with  an  absolute  average  deviation  of 
2.3%  for  63  data  points. 


Pressure  (bar) 

Figure  9  VE-P  projection  from  this  work  for  pseudo-binary  propane  in  66/34  wt% 
trans/cis- DHN  mixture  at  373  K  (filled  circles-65  mol%  C3,  filled  squares-55  mol%  C3,  filled 
triangles-45  mol%  C3,  unfilled  diamonds-35  mol%  C3,  unfilled  squares-25  mol%  C3,  unfilled 
triangles- 15  mol%  C3,  unfilled  circles-5  mol%  C3) 

Table  5  Tabulated  data  for  volume  expansion  of  C3  in  66/34  wt%  trans/cis- DHN  mixtures 


T|K) 

,  {mol  CS/md) 

VE,  P6S 

VE,  70  bar 

VE,  139  bar 

VE,  2Q8  aar 

323 

0.05 

1.014 

1.010 

1.012 

1.010 

0.15 

1.030 

1.032 

1.032 

1.023 

0.25 

1.050 

1.060 

1.064 

1.056 

0.35 

1.089 

1.084 

1.093 

1.080 

0.45 

1.131 

1.124 

1.126 

1.122 

0.55 

1.178 

1.169 

1.170 

1.173 

0.55 

1.277 

1.246 

1.235 

1.233 

373 

0.05 

1.007 

1.008 

1.011 

1.011 

0.15 

1.036 

1.037 

1.036 

1.035 

0.25 

1.080 

1.079 

1.076 

1.067 

0.35 

1.107 

1.110 

1.103 

1.09S 

0.45 

1.159 

1.156 

1.144 

1.137 

0.55 

1.212 

1.218 

1.210 

1.192 

0.55 

1.288 

1.306 

1.285 

1.262 

423 

0.05 

1.007 

1.017 

1.018 

1.007 

0.15 

1.047 

1.045 

1.041 

1.034 

0.25 

1.089 

1.081 

1.078 

1.074 

0.35 

1.125 

1.125 

1.117 

1.109 

0.45 

1.179 

1.193 

1.180 

1.165 

0.55 

1.278 

1.258 

1.248 

1.226 

0.55 

1.385 

1.386 

1.327 

1.314 

The  saturated  liquid 

curve  on  the  VE-P  isotherm  (solid  line)  shows 

an  increase  in  the 

volume  expansion  (VE)  with  pressure.  As  the  system  pressure  is  raised  for  a  two-phase  system, 

more  C3 

is  dissolved  in  the  liquid  resulting 

in  an  increase 

in  the  VE  of  the  liquid  phase.  The 

single-phase  VE  (dashed  curve)  is  relatively  insensitive  to  an  increase  in  pressure  over  the  range 
of  pressures  measured.  Even  though  the  mixture  density  decreases  with  increasing  C3 
concentration  (i.e.  the  fluid  is  less  liquid-like),  the  single  phase  mixture  is  relatively 
incompressible.  This  phenomenon  is  advantageous  at  the  high  pressures  required  to  obtain 
moderately  high  polymer  concentrations.  At  these  conditions,  the  expanded  fluid  phase  will 
exhibit  low  viscosities  and  enhanced  rates  of  mass  transfer  without  sacrificing  the  high  polymer 
concentrations  needed  to  obtain  a  highly  productive  hydrogenation  process. 

5.2  Viscosity 


Zeberg-Mikkelsen  and  coworkers  reported  the  viscosities  of  pure  isomers  of  DHN  over  a 
wide  range  of  pressures  and  temperatures  [28].  Having  obtained  three  pure  component 
parameters  for  both  the  cis  and  the  tnmv-DHN,  the  authors  rationalized  the  differences  in  the 
parameter  sets  using  differences  in  the  molecular  configurations  of  the  isomers.  The  same 
structural  features  that  make  the  cis  isomer  a  higher  density  fluid  were  cited  as  the  dominating 
factor  leading  to  significantly  higher  values  of  the  free  volume  overlap  and  lower  values  of  the 
characteristic  length. 

The  differences  between  saturated  liquid  viscosities  of  the  pure  DHN  isomers  and  our 
FVM  predictions  for  viscosity  are  shown  in  Figure  10.  The  66/34  wt%  trans/cis- DHN  saturated 
liquid  mixture  viscosity  is  also  shown  for  comparison.  The  model  results  shown  in  Figure  10 
were  computed  using  the  parameter  values  reported  by  Zeberg-Mikkelsen  et  al.  The  agreement 
between  the  FVM  predictions  for  the  pure  isomers  and  the  isomeric  mixture  of  DHN  is  excellent 
with  an  AAD  of  2.3%,  consistent  with  the  AADs  reported  by  Zeberg-Mikkelsen. 


Temperature  (K) 


Figure  10  Measured  saturated  liquid  viscosity  of  pure  isomers  of  DHN  and  a  66/34 
trans/cis-DHN  mixture  from  this  work,  Curves  are  FYM  predictions  using  SLEOS  for  density 
predictions 

Table  6  Tabulated  data  for  saturated  liquid  viscosity  of  isomers  of  DHN  and  66/34  wt% 


trans/cis-DHN  mixture 


323 

1.28 

1.86 

324 

1.45 

333 

1.11 

1.57 

348 

1.02 

343 

0.87 

1.34 

373 

0.75 

353 

0.66 

1.16 

388 

0.57 

363 

0.76 

1.02 

423 

0.45 

373 

0.66 

0.81 

383 

0.62 

0.81 

363 

0.56 

0.73 

403 

0.52 

0.67 

413 

0.47 

0.61 

423 

0.44 

0.56 

Isobaric  C3/DHN  compressed  liquid  mixture  viscosity  data  at  139  bar  for  six  C3 
concentrations  at  five  temperatures  are  shown  in  Figure  11.  Similar  measurements  were  made  at 
70  bar,  208  bar,  and  at  the  bubble  point  pressure  of  the  mixture.  The  FVM  was  fit  to  the 
C3/DHN  mixture  viscosity  data  for  all  three  isobars  by  varying  the  viscosity  binary  interaction 
parameter,  kac3-DHN,  as  shown  in  Equation  12.  When  the  binary  interaction  parameter  was  not 
used,  the  model  captured  the  data  well  qualitatively,  but  the  A  AD  was  14.5%  with  a  bias  of 
13.5%.  When  a  best-fit  binary  interaction  parameter  (kc3-DHN  =  -0.199)  was  used,  a  considerably 
better  fit  to  the  data  was  obtained  (AAD  of  5.7%  and  bias  of  0.29%). 


MIX 


Temperature  (K) 


Figure  11  Measured  viscosity  of  C3  in  66/34  trans/cis- DHN  at  139  bar  from  this  work, 
Curves  are  FVM  predictions  with  SLEOS  density  predictions  [kc3-DHN  =  -0.199] 


Table  7  Tabulated  viscosity  data  for  C3  in  66/34  wt%  trans/cis- DHN  mixtures 


0  mol%  C3  5  mol%  C3  15mol%C3  25mol%C3  35mol%C3  45mol%C3 


ETcSI 

nrai 

DtSl 

BTewa 

DtSl 

KZ31 

Klfl wl 

70 

324 

1.60 

325 

1.38 

324 

1.05 

324 

0.85 

323 

0.67 

323 

0.50 

348 

1.14 

349 

0.97 

348 

0.78 

349 

0.64 

348 

0.51 

348 

0.40 

373 

0.83 

373 

0.73 

373 

0.58 

373 

0.47 

373 

0.40 

373 

0.32 

398 

0.63 

399 

0.57 

398 

0.48 

399 

0.35 

398 

0.33 

398 

0.26 

423 

0.49 

423 

0.46 

423 

0.35 

423 

0.28 

423 

0.25 

423 

0.22 

139 

323 

1.75 

324 

1.49 

324 

1.17 

324 

0.90 

323 

0.71 

323 

0.53 

348 

1.25 

349 

1.09 

348 

0.85 

349 

0.69 

348 

0.55 

348 

0.43 

373 

0.91 

373 

0.80 

373 

0.64 

373 

0.52 

374 

0.43 

373 

0.34 

398 

0.68 

399 

0.62 

398 

0.52 

399 

0.38 

398 

0.36 

398 

0.28 

423 

0.54 

423 

0.50 

423 

0.39 

423 

0.32 

423 

0.29 

423 

0.23 

208 

50 

1.85 

324 

1.61 

324 

1.24 

324 

0.96 

323 

0.75 

323 

0.58 

75 

1.33 

348 

1.19 

348 

0.91 

348 

0.74 

348 

0.58 

348 

0.46 

ICO 

0.99 

372 

0.86 

373 

0.69 

373 

0.58 

373 

0.47 

373 

0.37 

124 

0.75 

399 

0.66 

398 

0.56 

398 

0.42 

398 

0.39 

399 

0.31 

149 

0.59 

423 

0.54 

423 

0.45 

423 

0.36 

423 

0.30 

423 

0.25 

The  influence  of  C3  on  the  mixture  viscosity  can  be  evaluated  with  the  viscosity 
reduction  (VR)  as  defined  by  the  expression, 


VR  =  1- 


r!MIX(T,P) 

tjdhn(T,P) 


(15) 


One  VR-P  isotherm  at  373  K  is  shown  in  Figure  12.  The  liquid  mixture  viscosity  measurements 
taken  on  the  two  phase  system  exhibited  similar  behavior  to  the  YE-P  isotherms.  The  VR  for  the 
single-phase  concentration  isopleths  is  mostly  insensitive  to  an  increase  in  pressure  over  the 
entire  range  of  the  data,  as  was  similarly  observed  in  the  volume  expansion  behavior.  The 
highest  VR  occurs  at  the  highest  C3  concentration  measured  (45  mol%  C3)  and  the  highest 


temperature  (423  K). 


0.90 


Pressure  (bar) 

Figure  12  VR-P  projection  from  this  work  for  mixtures  of  C3  in  66/34  trans/cis- DHN  at  373 
K  (unfilled  diamonds-saturated  liquid,  unfilled  triangles-45  mol%  C3,  unfilled  circles-35mol% 
C3,  filled  squares-25  mol%  C3,  filled  triangles- 15  mol%  C3,  filled  circles-5  mol%  C3) 

5.3  Correlation  of  viscosity  reduction  to  volume  expansion 

The  magnitude  of  viscosity  reduction  of  DHN  by  C3  at  various  temperatures  and 
pressures  is  almost  entirely  dependent  on  the  volume  expansion  of  DHN  by  C3.  We’ve  plotted 
the  viscosity  reduction  (necessarily  a  normalized  viscosity  value)  versus  the  volume  expansion 
(specific  volume  of  the  mixture  relative  to  pure  DHN)  in  Figure  13.  The  specific  volume  of  the 
mixture  captures  variations  in  temperature,  pressure,  and  dilution  thus  Figure  13  correlates  the 
viscosity  reduction  to  the  complete  physical  state  of  the  fluid. 


The  collapse  of  the  data  onto  a  single  curve  is  incomplete,  suggesting  that  other  physical 
attributes  of  the  mixture  also  influence  the  viscosity  reduction.  However,  the  degree  of 
separation  of  the  data  from  a  single  curve  is  relatively  small  therefore  much  of  the  viscosity 
reduction  can  be  attributed  to  the  expansion  of  the  fluid  phase  by  C3. 


Volume  Expansion 

Figure  13  VR  vs.  VE  from  this  work  for  mixtures  of  C3  in  66/34  /ra/is/dv-DHN 

6.  Conclusions 

A  range  of  temperatures,  pressures  and  compositions  of  C3/DHN  mixtures  suitable  for 
solvation  of  polystyrene  has  been  investigated  by  means  of  phase  equilibria,  density,  and 
viscosity  measurements.  Phase  equilibria  and  density  data  of  the  pure  components  and  mixtures 
were  captured  quite  well  by  the  SLEOS.  The  SLEOS  calculations  were  used  to  feed  biphasic 


and  single  fluid  densities  to  the  FVM  for  calculating  pure  component  viscosities  and  viscosities 
of  the  mixtures.  The  calculations  lead  to  very  good  agreement  between  the  FVM  for  viscosity 
and  the  data  with  a  single  adjustable  parameter.  Even  though  a  binary  interaction  parameter  for 
viscosity  was  required  to  obtain  a  good  fit  to  the  data  the  calculations  can  still  provide  utility  for 
future  calculations  that  require  interpolated  data  that  is  both  complex  and  laborious. 

List  of  symbols 

B  characteristic  of  the  free  volume  overlap  (FVM) 
c  total  number  of  components  in  the  mixture 

k  fitted  mixture  parameter  accounting  for  specific  binary  interactions  between  components 
in  the  mixture  (SLEOS  and  subscript  a  indicating  for  FVM) 
l  fitted  mixture  parameter  roughly  accounting  for  free  volume  differences  between 
components  in  the  mixture  (SLEOS,  not  used  in  this  work) 

M  molecular  weight 

n  moles 

P  pressure 

r  number  of  lattice  sites  per  molecule  (SLEOS) 

R  universal  gas  constant 

T  temperature 

V  specific  volume 

x  mole  fraction 

z  compressibility  factor 


Greek  symbols 


a  characteristic  barrier  energy  a  molecule  must  overcome  for  self-diffusion  to  occur  (FVM) 
s  interaction  energy  per  mer  (SLEOS) 

<f>  close-packed  volume  fraction  in  the  mixture  at  the  incompressible  state  (SLEOS) 

<p  fugacity  coefficient 

r|  dynamic  viscosity 

p  mass  density 

v  SLEOS  hard-core  molecular  volume  (SLEOS) 

£  characteristic  molecular  length  parameter  (FVM) 

Subscripts/superscripts 

a  indicating  binary  interaction  parameter,  k ,  is  for  FVM 
i,  j  for  compound  i,  j,  etc. 

characteristic  variable 
~  reduced  variable 

MIX  mixture  variable 

DHN  76/24  wt%  ?ra«.v/c/v-decahydronaphthalene  variable 

C3  propane  variable 

Abreviations 

VE  volume  expansion 

VR  viscosity  reduction 

C3  propane 

DHN  decahydronaphthalene  (assumed  76/24  wt%  trans/cis  unless  otherwise  stated  by  t  or  c 
indicating  the  trans  or  cis  isomer,  respectively) 
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